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Abstract 

A variational principle for Lagrangian densities containing derivatives of real order is for- 
rj^ , mulated and the invariance of this principle is studied in two characteristic cases. Necessary 
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and sufficient conditions for an infinitesimal transformation group (basic Nother's identity) 
are obtained. These conditions extend the classical results, valid for integer order deriva- 
tives. A generalization of Nother's theorem leading to conservation laws for fractional Euler- 
Lagrangian equation is obtained as well. Results are illustrated by several concrete examples. 
Finally, an approximation of a fractional Euler-Lagrangian equation by a system of inte- 
ger order equations is used for the formulation of an approximated invariance condition and 
corresponding conservation laws. 
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1 Introduction 

There are two distinct approaches in the formulation of fractional differential equations of models 
in various branches of science, e.g. physics. In the first one, differential equations containing 
integer order derivatives arc modified by replacing one or more of them with fractional ones 
(derivatives of real order). In the second approach, which we will follow in the sequel, one starts 
with a variational formulation of a physical process in which one modifies the Lagrangian density 
by replacing integer order derivatives with fractional ones. Then the action integral in the sense 
of Hamilton is minimized and the governing equation of a physical process is obtained. Hence, 
one is faced with the following problem: find minima (or maxima) of a functional 

C[u}= [ L{t,u(t), a D?u)dt, < a < 1, (1) 

J A 

where a Dfu is the left Riemann-Liouville fractional derivative, under certain assumptions on a 
Lagrangian L, as well as on functions u among which minimizers are sought. After this, the 
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Euler-Lagrange equations are formed for the modified Lagrangian, which leads to equations that 
intrinsically characterize a physical process. This approach has a more sound physical basis (see 
e.g. [H [SJ EH EH OH])- It provides the possibility of formulating conservation laws via Nother's 
theorem. Conservation laws are very important in particle reaction physics for example, where they 
are often postulated as conservation principles |37j . The central point is the fact that equations 
of a model, i.e., Euler-Lagrange equations, give minimizers of a functional. As a rule, fractional 
differential equations obtained through Euler-Lagrange equations contain left and right fractional 
derivatives, which makes this approach more delicate. 

In the last years fractional calculus has become popular as a useful tool for solving problems 
from various fields [2TJ [41] (see also |3[n|[THl[Till3lllll3lig|5S|3g|3S]). 

The study of fractional variational problems also has a long history. F. Riewe |32[ I3"3"] inves- 
tigated nonconservative Lagrangian and Hamiltonian mechanics and for those cases formulated 
a version of the Euler-Lagrange equations. O. P. Agrawal continued the study of the fractional 
Euler-Lagrange equations [IJ [21 13] ; for general fractional variational problems involving Riemann- 
Liouville, Caputo and Riesz fractional derivatives. G. S. F. Frederico and D. F. M. Torres [13] 
studied invariance properties of fractional variational problems with Nother type theorems by in- 
troducing a new concept of fractional conserved quantity which is not constant in time, so the term 
conserved quantity is not clear. Conservation laws and Hamiltonian type equations for the frac- 
tional action principles have also been derived in 39 . There are nowadays numerous applications 
of fractional variational calculus, see e.g. [3 El IE IS1 [IHl [HI [121 [IS1 IM [23 [23 [MJ [3D] - 

We can summarize the novelties of our paper as follows. First, we derived infinitesimal criterion 
for a local one-parameter group of transformations to be a variational symmetry group for the 
fractional variational problem ([!} . In previous work |13[ 1141 115] this was done only in a special 
case. Moreover, we separately consider cases when the lower bound a in the left Riemann-Liouville 
fractional derivative is not transformed (Subsection 12. ip , and when a is transformed (Subsection 
I2.2[) . Both cases have their physical interpretation. In the case when fractional derivatives model 
memory effects the lower bound a in the definition of derivative should not be varied. However, if 
fractional derivatives model nonlocal interactions (e.g. in nonlocal elasticity) the lower bound has 
to be varied. Second, we give a Nother type theorem in terms of conservation laws as it is done 
in the classical theory. This approach preserves the essential property of conserved quantities 
to be constant in time. Third novelty is the approximation procedure which is given for the 
Euler-Lagrange equations and infinitesimal criterion, therefore also for Nother's theorem. By the 
use of additional assumptions we obtain that appropriate sequences of classical Euler-Lagrange 
equations, infinitesimal criteria and conservation laws converge in the sense of an appropriate space 
of generalized functions to the corresponding Euler-Lagrange equations, infinitesimal criterion and 
conservation laws of the fractional variational problem ([!} . There are also a number of illustrative 
examples which complete the theory. 

The paper is organized as follows. To the end of this introductory section we provide basic 
notions and definitions from the calculus with derivatives and integrals of any real order, and 
recall the Euler-Lagrange equations of ([T]) for (A, B) C (a, b) (so far only the case (A, B) = (a, b) 
was treated, cf. [U 12])- Section [2] is devoted to the study of variational symmetry groups and 
invariance conditions. In Section [3] we derive a version of Nother's theorem which establishes a 
relation between fractional variational symmetries and conservation laws for the Euler-Lagrange 
equation, which have the important property of being constant in time. We use in Section [4] the 
approximation of the Riemann-Liouville fractional derivative by a finite sum of classical derivatives. 
In this way and by introducing a suitable space of analytic functions we obtain approximated 
Euler-Lagrange equations, approximated local Lie group actions with corresponding infinitesimal 
criteria and Nother theorems. 

1.1 Preliminaries 

In the sequel we briefly recall some basic facts from the fractional calculus. 

For u € ^([a, b]), a > and t £ [a, b] we define the left, resp. right Riemann-Liouville fractional 
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integral of order a, as 

J?u=-^- { (t-0) a - l u{0)d6, resp. t I£u = / (0 - t)"" 1 ^) t». 

r ( a ) 7 a r (a) 7t 

Left, resp. right Riemann-Liouville fractional derivative of order a, < a < 1, is well defined for 
an absolutely continuous function u in [a, 6], i.e. u € AC([a, 6]), and i G [a, b] as 

aD t u--J t U -^—- ) -J ajr - W ;d6, (2) 

resp. 

t D?u=(-4-) al-u = ( 4-) f de. 

b \ dt) b r(l-a)V dt) Jt (0-t) a 

If /, g G AC ([a, b]) and < a < 1 the following fractional integration by parts formula holds: 

[ f(t) a D?gdt= f g(t) t D?fdt. (3) 

J a J a 

For the left Riemann-Liouville fractional derivative © of order a (0 < a < 1) we have: 

a ZX Q w= — r/ ? ^ rf6* H K ' tea, , 4) 

4 r(l - a) 7 a (t - 0) Q r(l-a)(i-a) a 1 J w 

The integral on the right hand side is called the left Caputo fractional derivative of order a and 
is denoted by c a Dfu. Similarly, the right Caputo fractional derivative t^b u ^ s defined. 

It follows from {4| that the left Riemann-Liouville equals to the left Caputo fractional derivative 
in the case u(a) = (the analogue holds for the right derivatives under the assumption u(b) = 0). 
The same condition, i.e. u(a) = 0, provides that -^ a D"u = a D^^u. 

We introduce the notation which will be used throughout this paper: derivatives of Lagrangian 

dL 

L = L(t,u(t), a Dfu) with respect to the first, second and third variable will be denoted by 

■7— and — — ~~~ — , or by d\L, d^L and d^L respectively. 
ou o a D"u 

1.2 Euler-Lagrange equations 

As stated in Introduction, we solve a fractional variational problem 

C[u]= [ L(t,u(t), a D?u)dt -> min, < a < 1. (5) 

J A 

(A, B) is a subinterval of (a, b), and Lagrangian L is a function in (a,b) xlxK such that 

L eC x ((a, b) xRxM) ] 

and \ (6) 

t ^ d 3 L(t,u(t), a D?u) e AC([a,b\), for every u G AC([a,b}) J 

Solutions to (JS|) are sought among all absolutely continuous functions in [a, b], which in addition 
satisfy condition u(a) = ao, for a fixed ao G R. 

One can consider more general problems with Lagrangians depending also on the right Riemann- 
Liouville fractional derivative. Our results can be easily formulated in that case, and because of 
simplicity we shall consider only Lagrangians of the form <j5j> - Moreover, in ([5]) one can consider 
the case a > 1, but this also does not give any essential novelty and because of that it is skipped. 

We will recall results about the Euler-Lagrange equations obtained in [TJ [2J H] . As we said, we 
consider the fractional variational problem defined by ([5]) . 
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Fractional Euler-Lagrange equations, which provide a necessary condition for extremals of a 
fractional variational problem, have been recently studied in [TJ [H 0]- Let A = a and B = b. 
Euler-Lagrange equations are obtained in [J [5] : 

dL / dL 



In terms of the Caputo fractional derivative, the Euler-Lagrange equation ([7]) is given in [4] and 
reads: 

dL 



c n a 



( dL ) 


dL 


1 1 


\d a D?u) 


+ d a D?u 


t=b T(l-a) (b-t) a 



du 

Euler-Lagrange equations for ([5J are derived in [4] : 



= 0. (») 



MiS>)--<im-J - Ms(o " 4) - (10) 

Equation @ is equivalent to 

Remark 1. If we replace the Lagrangian in ([5]) by L(t,u(t), a D"u,u(t)) then Euler-Lagrange 

dL 

equations will be also equipped with a 'classical' term — on the left hand side. It can be 

at dii 

shown that in this case infinitesimal criterion and conservation law can be obtained by combining 
the results of the present analysis (see Theorems El [11] and [15]) and the classical theory (see e.g. 

2 Infinitesimal invariance 

Let G be a local one-parameter group of transformations acting on a space of independent and 
dependent variables as follows: (t, u) = ■ {t,u) — (H^(i, it), ^(t, u)), for smooth functions 
and ^f v , and g v £ G. Let 

v = T[t,u)— +£{t,u) — 
dt du 

be the infinitesimal generator of G. Then we also have 

t = t + r]T(t,u)+o(r]) ^ 

U = U + 7]^(t, U) + 0(77). 

We introduce the following notation (cf., e.g., [3D]): At = ^|„ =0 (i—t) and Am = ^-|„ =0 (u(t) — 
m(tj)). More precisely, the notations At and Am denote lim^o *^^~ t and lim^o w f^~ wl ^> ; 
respectively. It follows from (|12p that A< = r, and writing the Taylor expansion of u(t) = 
u(i + r\T(t, u) + 0(77)) at 77 = yields that Am = £. 

On the other hand, if we write Taylor expansion of u(t) at t = t we obtain 



d 

Am = — 

or] 



(u(t) - u{t)) + uAt, 



?7=0 



or, if we introduce the Lagrangian variation 

r d 

du := — 

drj 



(u(t)-u(t)) (13) 

?? = 



4 



we obtain 
Thus, we have 



Au = Su + uAt. 
5u = £ — TU. 



Note that St = 0. 

In the same way we can define AF and SF of an arbitrary absolutely continuous function 
F = F(t,u(t),u(t)): 



AF 



d 
drj 

5F = 



(F(t,u(t),u{tj) - F{t,u(t),u(t)) 

77=0 ^ 



dF A 



dF . OF A . 
— Au + — Au 

au au 



dt] 



/ . \ OF dF 

F{t, u(t), ti t - F(t, ti t), u(t = — Su + —Sii, 
, V ) au au 



and 



AF = SF + FAt. 



It will be apparent in the forthcoming sections that, if we want to find an infinitesimal criterion, 
we need to know A a Dfu and AC. Therefore, we have to transform the left Riemann-Liouville 
fractional derivative of u under the action of a local one-parameter group of transformations (|12j) . 
There are two different cases which will be considered separately. In the first one, the lower bound 
a in a -Df u is not transformed, while in the second case a is transformed in the same way as the 
independent variable t. Physically, the first case is important when Q D"u represents memory 
effects, and the second one is important when action on a distance is involved. 



2.1 The case when a in a Dfu is not transformed 

In this section we consider a local group of transformations G which transforms t G (A, B) into 
i S (A, B) so that both intervals remain subintervals of (a, 6), but the action of G has no effect on 
the lower bound a in a D"u, i.e., r(a, u(a)) = 0. So, suppose that G acts on i,u and a D"u in the 
following way: 

g v ■ (t, u, a D?u) := (i, u, a Df u), 
where t and u are defined by (|12|) . In this case we have: 



Lemma 2. Let u G AC ([a, b]) and let G be a local one-parameter group of transformations given 
by (E§. Then 

A a D?u = a D?Su + i a D?u ■ r(t, u(t)), 



dr 



where 



S a D?u 



d_ 

drj „ = 



?7 = >> 



u- a Dfu) = a D?5u. 



Proof. To prove that 5 a Dfu = a DfSu it is enough to apply the definition of the Lagrangian 
variation (|13[) . Also by definition we have 



A a D?u 



Thus, 

A a D?u = 



1 



T(l - a) dr] „ =0 



d 

drj n= 
1 u{6) 



a Dfu - a D?u 



= a D«8u+ — 



a D?5u 



dr] v=0 
d di 



dij a (i-9) 
a Dfu - a D?u 

( a Dfu - a-D™ 



d6± — 



dtj a (t-ey 



■d6-± 



u(0) 
dtj a (t-6)° 



dt dr] 

= a D?5u+j ta D?u-T(t,u(t)) 
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□ 



Lemma 3. Let C[u] be a functional of the form C[u] = J A L(t,u(t), a Dfu) dt, where u is an 
absolutely continuous function in [a,b], (A, B) C (a,b) and L satisfies {6p. Let G be a local 
one-parameter group of transformations given by I112\) . Then 



where 5C = f b SLdt. 

J a 



AC = 5C + (LAt) 



Proof. Again it is clear that SC = 5L dt. For AC we have 



drj 



L(t,u(t), a D?u)dt 



L{t,u{t), a D?u)dt 



d_ 

drj „ =0 
d_ 

drj „ =0 



L(t, u(t), a Dfu) (1 + 777-0, u(t))) dt 



L(t,u(t), a D?u)dt 



L{t,u{f), a Dfu)dt 



L(t, u(t), a Dfu) f (t, u{t)) dt. 



L(t,u(t), a D?u)dt 



This further yields 
AC : 



AL(t,u(t), a D?u)dt + / L{t,u(t), a D?u)i(t,u(t))dt 



8L{t,u{t), a D?u)dt- 



dt 



Lit, u(t), a D?u) r(t, u(t)) dt 



L(t, u(t), a D?u) f(t, u(t)) dt 
/ SLdt+(Lr) 

J A A 



and the claim is proved. 



□ 



We define a variational symmetry group of the fractional variational problem (JSJ, which we 
call a fractional variational symmetry group: 

Definition 4. A local one-parameter group of transformations G Alty) is a variational symmetry 
group of the fractional variational problem |5p if the following conditions holds: for every [A' , B'] C 
(A, B), u = u(t) £ AC ([A', B']) and g v e G such that u(t) = g v ■ u{t) is in AC([A',B']), we have 



L(t,u(t), a Dfu)dt 



B' 



L(t,u(t), a D?u)dt. 



(14) 



We are now able to prove the following infinitesimal criterion: 



Theorem 5. Let C[u] be a fractional variational problem f5]) and let G be a local one-parameter 
transformation group with the infinitesimal generator v = r{t,u)dt + £(t,u)d u - Then G is a 
variational symmetry group of C if and only if 



dL 



d a D?u 



Lf = 0. 



(15) 



G 



Proof. Suppose that G is a variational symmetry group of C. Then (TlH) holds for all subintervals 
(A',B') of (A, B) with closure [A',B'] c(A,B).We have: 



AC = [ 5Ldt + (LAt) 

J A' 



A 
B' 

A' 
B' 

A' 
B' 



dL . dL 



!Tu 5u + lWu 5{aD>)],U + [LX ' ] 

— {Au - uAt) + - — a D?(Au - uAt) ) ,11 -:- [LAI ) 
B ' / fit bt F)T 



Applying the Leibnitz rule for f t {LAt) we replace ^At+^uAt by |(LAf)- 3 ^„ M (^ a Dfu)At- 
LAt^h Then 



dt du d a Dfu 



<>L "^ a D?u)At - a D?(uAt)) + LAt {1) ) dt. 



d a D?u\\dt 



Since AC has to be zero in all (A' , B') with [A 1 , B'\ C (A, B), the above integrand has also to be 
equal zero, i.e., 



where we have used that Au = £ and At = r. Hence, necessity of the statement is proved. To prove 
that condition (p~5|) is also sufficient, we first realize that if (fT5j) holds on every C (A,B), 

then A£ = in every [A , £?'] C (A, £?). Thus, integrating AC from to ij we obtain (fl4|) . for rj 
near the identity. The proof is now complete. □ 

In the following example we calculate the transformation of the fractional derivative a Dfu 
under the group of translations. In Subsection 12.21 we shall perform such calculation also in the 
case when a in a D^u is transformed. 

Example 6. Let G be a local one-parameter translation group: (i, u) = (t + r),u + rf), with the 
infinitesimal generator v = d t + d u . Then u(t) = u(t — rf) + r\ and by a straightforward calculation 
it can be shown that 

2.2 The case when a in a Dfu is transformed 

Now let the action of a local one-parameter group of transformations (|T2|) be of the form 

g n ■ (t,u, a D"u) := (f,n, g £)fu). 

Thus, the one-parameter group also acts on a and transforms it to a, where a = g v ■ t\ t=a - This 
will also influence the calculation of A a Dfu and AC: 
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Lemma 7. Let u G AC([a, b]) and let G be a local one-parameter group of transformations flty) . 
Then 



A a D?u = a D?Su + j e D a t u ■ r(t, u(t)) + — (f " ( ° ) a+1 r(a, «(o)), 



(16) 



where 



8 a D?u = — 

drj „ =0 



Proof. Again using (|T3l) we check that 8 a Dfu — a Df8u. To prove (|16|) we start with the 



definition of A a D? 



Thus, 



A Q D t Q u = — 
dry 



8 Z>f u - a Dfu 



A a D?u 



1 



T(l - a) dr) „ =0 
1 



u(0) 



diJa (* 



■d<9± 



d(9 



d*A (t-0) a di7 a (t 



u(0) 



-de 



a Df8u 



T(l - a) d?? ,, =0 



d r u(6») 



■<20 + 



u(0) 



d0 



d(9. 



Note that in the last term we can interchange the order of integration and differentiation with 
respect to t. This gives 



a m 

r(l- a) dr] , 1=0 J- a (i-6) a 



■d9. 



If we differentiate this integral with respect to rj at r\ — we eventually obtain (fTBj) . 



□ 



Lemma 8. Let C[u] be a functional of the form C[u] = jf L(t,u(t), a Dfu) dt, where u is an 
absolutely continuous function in (a,b), (A, B) C (a, b) and L satisfies {TjJ). Let G be a local 
one-parameter group of transformations given by Then 



A£ = 5£+ (LAt) 



<9L 



daDi 



-dt. 



where SC — SL dt. 

Proof. Clearly, SC — f SL dt. On the other hand we have 



dn _ = 



( [ L(i,u(t),- a D?u)dt- [ L(t,u(t), a D?u)dt) 

.7=0 \ JA J A / 



d 

dr) 



i=o \ J A 



L(t,u(t),aDfu)(l + rjf(t,u(t)))dt- / L(t,u(t), a Dfu) dt 



ALdt+ / L(t,u(t), a D"u)f(t,u(t))dt 



(17) 



/dL dL dL \ f B 

(—At+—Au+——A a D?u)dt + / L(t,u(t), a D?u)i(t,u(t))dt. 
V at au o a D"u / J A 
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We now apply (IT5|) to obtain 



AC = ^(^At+^u + ^ + ^^+A^.^^)) 



+ r(l- a ) (t-a)°+i T(a ' U{a)) ) ) M + I A L{t ' aD>) ^ U(t)) dt 



B dL , 
dt, 



i a d a D?u 

f£ + H, + r(l-a)(t-a)^ ^ tt ^A ^> 



dt. 



□ 



Remark 9. It should be emphasized that the last summand on the right-hand side of (|16p. as 
well as of (|17p. appears only as a consequence of the fact that the action of a transformation 
group affects also a, on which the left Riemann-Liouville fractional derivative depends. In the 
case when u(a) = or a — 1 that term is equal to zero (since lim,-^!- T(l — a) = oo) and 

AC = SC + (LAt) . 

A 

Now we define a variational symmetry group of the fractional variational problem ([5]) as follows: 

Definition 10. A local one-parameter group of transformations G is a variational symmetry 
group of the fractional variational problem {5p if the following conditions holds: for every [A',B'] C 
(A, B), u = u(t) £ AC([A' , B']) and g v e G such that u(t) = g v ■ u(t) is in AC ([A 1 , B']), we have 

pB' r^' 
/ L(i,u(t), a D?u)di= L(t,u(t), a D?u)dt. 

J A' J A' 

Recall that we are solving the fractional variational problem ([5]) among all absolutely continuous 
functions in [a, b], which additionally satisfy u(a) — 0. Therefore, the last term in both (|16[) and 
(|17| vanishes, and the infinitesimal criterion reads the same as in Theorem [5j 

Theorem 11. Let C[u] be a fractional variational problem defined by {5|] and let G be a local one- 
parameter transformation group defined by $Wj) with the infinitesimal generator v = r(t,u)dt + 
£(t, u)d u . Assume that u(a) = 0, for all admissible functions u. Then G is a variational symmetry 
group of C if and only if 

Proof. See the proof of Theorem [5] □ 

Example 12. Let v = d t + d u be the infinitesimal generator of the translation group (t, u) = 
(t + 7/, u + rf). Then u(t) = u(t — rf) +rj and 

- a D?u = a D?{u + rj). 

3 Nother's theorem 

In the formulation of Nother's theorem we shall need a generalization of the fractional integration 
by parts. 
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Lemma 13. Let f. g £ AC ([a, b]). Then, for all t G [a, b] the following formula holds: 



f(s)- s DZgds 



J a D a J-g{s)ds 

9(b) g(t) 



r(i-a) 



(19) 



_(fe-s) a (*-s) Q 7* (o--s) a 

Proof. In order to derive (fT9")) we shall use the representation of the right (resp. left) Riemann- 
Liouville fractional derivative via the right (resp. left) Caputo fractional derivative (QJ, as well as 
©: 



f{s)- s D%gds = 



i 



/(«) 



r(l-a) 
1 

' r(l-a) 

6 gir) 



9(b) 
(b-s) 

9(b) 



(<? 



(a- s y 



da± 



{b - s) a 
9(t) 



9{°) 



da 



■ da 



ds 



(t 



ds 



f(s) ■ s D?g + f(s) 



T(l - a) 



9(b) 
(b-sY 



g(t) 

(t - s) a 



(a - s) a 



da 



ds 



a D°f-g(s)+f(s) 



r(i - a) 



9(b) 
(b - s) a 



9(t) 

(t - s y 



9±) 
(a - s) a 



da 



ds. 



□ 



Remark 14. If we put t = b in (fT9"|) . we obtain Q. 

The main goal of symmetry group analysis in the calculus of variations are the first integrals 
of Euler-Lagrange equations of a variational problem, that is a Nother type result. 

Analogously to the classical case, an expression 4:P[t, u(t), a Dfu) = is called a fractional first 
integral (or a fractional conservation law) for a fractional differential equation F(t, u(t), a D"u) = 0, 
if it vanishes along all solutions u(t) of F. 

In the statement which is to follow, we prove a version of the fractional Nother theorem. As we 
will show, a fractional conserved quantity will also contain integral terms, which is unavoidable, 
due to the presence of fractional derivatives. If a = 1 then (|20[) reduces to the well-known classical 
conservation laws for a first-order variational problem. 

Theorem 15. (Nother's theorem) Let G be a local fractional variational symmetry group defined 
by of the fractional variational problem and let v = r(t,u(t))dt + £,(t,u(t))d u be the 
infinitesimal generator of G. Then 

Lr + jX^-^^-^-ur) s D % (^))ds= const. (20) 

is a fractional first integral (or fractional conservation law) for the Euler-Lagrange equation (0). 
The fractional conservation law \20) can equivalently be written in the form 



Lt 



1 



T(l - a) 
_d_ 

da 



8 3 L(B, u{B), a D%u) d 3 L{t, u(t), a D?u) 



(B - s) a 



b ±d 3 L(a,u(a), a D«u) 



(a - s) a 



da 



(t- s y 

ds = const. 



(21) 
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Remark 16. This theorem is valid in the case when the lower bound a in a Dfu is not transformed, 
as well as in the case when a is transformed. In the second case we in addition have to suppose 
that u(a) — 0, for all admissible functions u, which provides the same form of the infinitesimal 
criterion in both cases. 



Proof. We want to insert the Euler-Lagrange equation (JTTJ) into the infinitesimal criterion ([15 
or (fT8| . Hence, we will write the latter in a suitable form: 



dL dL . dL d \ .dL ^ Q .^ . . dL 



d 



dt 



The middle term in the last equation is the Euler-Lagrange equation (fTTj) multiplied by £ — iir. 
Therefore, 

|(ir + jf (jj? « - - ({ - -).^(^>) *) - o 

for all solutions of the Euler-Lagrange equation, and (f2T))) holds. If we now apply Lemma IT51 we 
can transform {20) into ([H]). □ 

In the next example we consider the case in which the Lagrangian of a fractional variational 
problem does not depend on t explicitly. 

Example 17. Consider a variational problem whose Lagrangian does not depend on t explicitly: 

C[u] = [ L{u(t), a D?u)dt. (22) 



As in the classical case, the one-parameter group of translations of time is a variational symmetry 
group of (|2"2"1) . Indeed, t = t + r],u = u (hence u(t) — u(t — rj)) and g,Df u = a Dfu (see Example 
□2J). Therefore, 

rB rB 

L(u(t),aDfu)di= / L(u(t), a D?u)dt. 



This fact can be also conhrmed by the infinitesimal criterion: 

dL 

where the last equality holds since u(a) = 0. Thus, j^ a Dfu = a D"-^u (see Introduction). 

Using Nother's theorem we may write a fractional conservation law which comes from this 
translation group: 

l + 1 ( - ^ -i&f-u +ii - tD% (^) ) ds = const - 



or equivalently 



t 

ii(s) 



I \d 3 L(u(B), a D%u) d z L{u{t), a D?u) 



r(l-a) 



» £d 3 L(u(a), a D«u) 
(a - s) a 



ds = const. (23) 
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Note that (gSJ) tends to 



dL(u(t),u(t)) 

L — u(t) • = const. 

ou 



when a — > 1 , which is the classical energy integral (see [40]V 

Example 18. As a concrete example related to Example [171 let us consider a fractional oscillator 
for which the Lagrangian reads: 



1 v 2 
L=-( a D?u) 3 -«?-, (24) 

where u is a constant (frequency). The fractional variational problem consists of minimizing the 
functional Ldt, among all smooth functions which satisfy the initial conditions u(0) = and 
u^'(0) = 1. The Euler-Lagrange equation for such an L is: 

uj 2 u - t D? D?u = 0. (25) 

Since L does not depend on t explicitly and u(0) = 0, the translation group (t, u) — (t + r), u) is a 
fractional variational symmetry group for (|24p. It generates the following fractional conservation 
law for (|2"5|) : 

i 2 ft 

-( a D?u) 2 -w 2 y + J ( - D>« ■ Dy + • s £>?( £>») = const. 

4 Approximations 

In this section we use approximations of the Riemann-Liouville fractional derivatives by finite 
sums of derivatives of integer order, which leads to variational problems involving only classical 
derivatives. We examine the relation between the Euler-Lagrange equations, infinitesimal criterion 
and Nother's theorem obtained in the process of approximation and the fractional Euler-Lagrange 
equations, infinitesimal criterion and Nother's theorem derived in previous sections. 
In the sequel we make the following assumptions: 

(i) L G C N ([a, b] x R x R), at least, for some N G N. 

(ii) We will simplify the following calculations by considering the case A = a and B = b. 

(iii) Let (c, d), —oo < c < d < +oo, be an open interval in R which contains [a, b], such that for 
each t G [a, b] the closed ball B(t, b — a), with center at t and radius b — a, lies in (c, d). 

In addition we shall separately consider two cases: 

(a) Let u G C°°([a, b]) such that u(a) — ao> u(b) — b , for fixed ao,6o G R, and L 3 (where I/3 
stands for d^L) be a function in [a,b] defined by t h-> £ 3 (i) = L 3 (t,u(t), a Dfu), t G [a, 6]. 
Let L3 (6, 60, = 0, for all i G N, meaning that for (t,s,p) 1— > L 3 (t,s,p), t G [a, &], s,p G R, 
the following holds: 

-^±{b,b ,p)=Q; -^-l(6,6 ,p)=0; -^(6, b ,p) = 0, Vp G R. 

(b) Let u G C°°([a, &]) such that w w (6) = 0, for all i 6 No, and it(a) = a , for fixed a G R. 
Let L^\b) = £3^(6,0, aDftii) — 0, for all i G N and for every fixed u, meaning that for 
(t, s,p) H> L 3 (t, s,p), t G [a, 0], s,p G R, the following holds: 

^(&,0,p) = 0; ^(6,0,p) = 0, VpGt. 
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Let / be a real analytic function in (c,d). Then according to 36, (15.4) and (1-48)] we have 
that 

'a\ (t — a)' 1 



aD"f = FFZi \f®®> teB(t,b-a)c(c,d), (26) 

rzk W r(* + l-a) 



cA (-l^aT^ - a) 



where , , 

r(l-a)r(i + l) 

Consider again the fractional variational problem ([5]) . Assume that we are looking for a mini- 
mizer u G C 2Ar ([a, 6]), for some iV G N. We replace in the Lagrangian the left Riemann-Liouville 
fractional derivative a Dfu by the finite sum of integer- valued derivatives as in (|26[) : 

b L(t, u(t),J2 (") r(* + i-l) " w W) ) dt = jf Z(< ' «(*).« w (*). « (a) (*),..-, *• 

(27) 

Now the Lagrangian L depends on t, u and all (classical) derivatives of it up to order N. Moreover, 

d 3 L, . . . , d N+2 L G C N ~ l ([a, b] x R x K), since d.Z = d 3 £(°) r(7+f-J) » 4 = 3 ' ^ + 2- For the P rob " 
lem (|27|> we can calculate the Euler-Lagrange equations, infinitesimal criterion and conservation 
laws (via Nother's theorem). 

Now we will consider the fractional variational problem (JS|) in the case (a) and in the case (b) . 

The main result of this section is stated in the following theorem. 

Theorem 19. Let C[u] be a fractional variational problem 0) which is being solved in the case (a) 
or (b). Denote by CL the fractional conservation law &20\), and by CLn the fractional conservation 
law for the Euler-Lagrange equation t28\) which corresponds to the variational problem Then 



CLjsi — » CL in the weak sense, as N — » oo. 

Convergence in the weak sense here means that as a test function space we use the space of 
real analytic functions as follows. (Note that analytic functions are chosen in order to provide 
convergence of the sum as in (|26|) .) 

Let A((c, d)) be the space of real analytic functions in (c, d) with the family of seminorms 

P[m,n]i<P) ■= Sup \<p(t)\, !peA((c,(l)), 

[m,n] 

where [to, n] are subintervals of (c, d). Every function / G C([a, 6]), which we extend to be zero in 
(c, d)\[a, b], defines an element of the dual A'((c, d)) via 

<P»(f,<p)= [ f(t)<f(t)dt, <peA((c,d)). 

J a 

Before we prove the main Theorem 1191 we shall prove several auxiliary results. 

First we recall a result from [4] which provides an expression for the right Riemann-Liouville 
fractional derivative in terms of the lower bound a, which figures in the left Riemann-Liouville 
fractional derivative. 

Proposition 20. Let F G C°°([a,b]), such that F^(b) = 0, for all i G N 0j and F = in 

(c, d)\[a, b]. Let tD^F be extended by zero in (c, d)\[a, b]. Then: 

(i) For every i G N, the (i — l)-th derivative of t h- > F(t)(t — a) l ~ a is continuous at t — a and 
t = b and the i-th derivative of this function, for i £ No, is integrable in (c, d) and supported 
in [a, b] . 

(ii) The partial sums Sn, N G No, 

dyf I a\ (t - a) 



I 0, t G (c, d)\[a, b], 

are integrable functions in (c, d) supported in [a, 6]; 
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(Hi) 



(-- 

\ dt 



tD%F = ^l --) ( F ■ 

i=0 



a\ (t-a) 1 



i J T(i + 1 - a) 



in the weak sense. 

Proof. See H Prop. 4.2]. □ 
The Euler-Lagrange equation for (|2"T|) is of the following form: 

N 



^ ( dt) On" 

4=0 



idL =0. 



This is equivalent to 



Remark 21. The Euler-Lagrange equation (|28p provides a necessary condition when one solves 
the variational problem (|27l) in the class C 2N ([a, b]), with the prescribed boundary conditions at 
a and b, i.e., u(a) = oq and = bo, oq, &o & re fixed real numbers. 



The following theorem shows that the Euler-Lagrange equation (|28|) converges to ([8]), as N — > 
+oo, in the weak sense. To shorten the notation, we introduce Pjv and P for the Euler-Lagrange 
equations in (|28|) and ((8|) respectively. 

Theorem 22. Let C[u] be the fractional variational problem $5$ to be solved in the case (a) or 
(b). Denote by P the fractional Euler-Lagrange equations (G|) ; and by Pn the Euler-Lagrange 
equations i28\) , which correspond to the variational problem {27}) , in which the left Riemann- 
Liouville fractional derivative is approximated according to H26}) by a finite sum. Then in both 
cases a) and b) 

Pn — >• P in the weak sense, as N — > 0. 
Proof. See % Th. 4.3]. □ 
Let 

d d 

vw = t n — +£ntt (29) 
ot ou 

be the infinitesimal generator of a local one-parameter variational symmetry group of (|27[) . The 
Af-th prolongation of is given by 

N „ 
(JV) . ~ O 



with 

d l 



N -TNUW) +T N U<- 1+1 \ 



Remark 23. It should be emphasized that the vector field v^v in (|2"9")l differs (in general) from 
the one introduced at the beginning of Section [2j 

Theorem 24. Let C[u] and C[u] be fractional variational problems f5j) and \27ty respectively. 
Denote by LC and ICm the corresponding infinitesimal criteria. If tn — > t and £jv £ as 
N — > oo, uniformly on compact sets, then 

ICn —> LC uniformly on compact sets, as N — > oo. 
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Proof. The infinitesimal criterion for the variational problem (|27[) says that a vector field vjv 
generates a local one-parameter variational symmetry group of the functional (|27[) if and only if 



prWv N (L)+Lt N = 0. 



This is equivalent to 



dL dL 
= TN-rr- + S.N Tr- 
ot Ou 



N / ji 



if 



AT 



N , 



i=0 



cA (z — a)(t — a) 



i — a — 1 



r(i + 1 - a) 



a\ (t — a)' 



oj r(i - a) 



AT 



«) +rA r ^ +1 M^ 



i=l 

tnOiL + £,n02Lv 

N 



a\ (t- a y- a 



E 



(Cat -tnu 



a^w , M (*-o)- a 
j + Wr(i-a) 



&i 



± 



TATM^ 1 ' + T N 



OJ r(i-a) 



N 

E 

i=0 



oA, (i-a)(t-o)* - ° 
+ 1 - a) 



AT 

E 



a\ (i — a) 1 



J V»/r(i + l-a) 

A 7 " 



c^L + if at 



T N d\L + ^at9 2 L 



aA (t-a)^ Q 



i/ T(i + 1 - a) 



A ; 



aA (i - a) l ~ a 



i J T(i + 1 — a) 



i TAr 



<9 3 L + it 



A 



(30) 



So, if tat — >• r and £ at — >■ £ as — ► 00, then the last expression on the right hand side tends to 

rd x L + id 2 L + (aD?{t - ut) + (^ a £) t a u)r)fti + if, 
and that is exactly the infinitesimal criterion (|15|) for §5$). □ 



Proof of Theorem 1191 First, we will derive the form of a conservation law for the Euler- 
Lagrange equation (|28|) of the approximated variational problem (|27[). which comes from a varia- 
tional symmetry group with infinitesimal generator (|29|) . For that purpose we need to insert the 
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Euler-Lagrange equations (1251) into the infinitesimal criterion (J3UJ): 



N 

E 



a\ (t-a) 1 



i J T(i + 1 — a 



dt 



N 



i=0 



a\ {t~ a y- a 



iJT(i+l-a) 



d 3 L + Ltn ± TjfU^foL 



TN 



N 



a\ (t-a) l - a 



N / x 



if. 



TV 



i=0 

d 



a\ (t-a) l - a 



±(& _ Trt a» )E (_| 



iy r(i + 1 - a) 

a\ (t — a) 1 



i=0 



i y r(i + 1 — a) 



JY 



dy ( f a\ (t-a) 



dt 



i)T(i+l-a) 



dt 



N 



i=0 



i J T(i + 1 — a) 



N , 

■to-W')E(-s)' 



dy((a\ (s-a) l ' a 



i J T(i + 1 — a) 



ft»L 



We recognize that the expression in brackets in the middle term in this sum is the Euler-Lagrange 
equation (|28|) . and hence vanishes. Therefore, we obtain that the following quantity is conserved: 



N 

E 



a\ (s- a y- a 

i Jr(j + l-a) 

N 



i=0 



cA (s — a) 1 



iy r(i + 1 - a) 



(31) 



ds = const. 



This of course is one way to write the first integral of the Euler-Lagrange equations (|28|). which 
corresponds to the local one-parameter variational symmetry group generated by (|29|) . 

Applying Proposition [20] (ii) we obtain that (|3Tj) converges to ([20]) in the weak sense, provided 
that {tn 7 £,n) —> ( T ,£)> as N —> oo. □ 



Remark 25. Numerical analysis of fractional differential equations by the use of Theorems [191 1221 
and[M]is not established and here we pose an open problem related to the use of our approximation 
procedure in applications. 



5 Concluding remarks 

In this paper we studied variational principle containing fractional order derivatives. Especially, 
we obtained a condition which the symmetry group for this principle has to satisfy. This symmetry 
group is used for the formulation of a conservation law which generalizes the classical Nother's 
theorem. In our work we used left Riemann-Liouville derivatives in Lagrangian density only. An 
extension including both left and right fractional derivatives could be easily performed. Using 
a method proposed in our previous work (see [1]) we approximated fractional Euler-Lagrangian 
equation with finite systems of integer order differential equations. In this way we obtained invari- 
ance conditions and corresponding conservation laws for such systems so that these conservation 
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laws converge to the fractional conservation laws when the number of integer order equations tends 
to infinity. 
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